RANDOM ANTENNA ARRAY INTERFEROMETER FOR RADIO LOCATION 
FIELD OF THE INVENTION 

This invention relates to a method and system for the accurate localization of multiple 
radio frequency signals over a large area. The innovation presented here was 
5 developed in the context of miniature (<1gram) wildlife tracking devices, but is 
generalized to the problem of locating either CMDA or non-CDMA transmitters. 

BACKGROUND 

Despite the revolutions represented by GPS technology and advanced direction-of- 
arrival algorithms (DoA), the problem of determining precise (< 1m resolution) 

10 positions of wildlife radio transmitters remains to date essentially unsolved. There 
are a number of reasons for this, principal among them being the low-duty cycle 
nature of the signals (required by practical energy/payload constraints), which 
prevents coherent demodulation. Payload constraints (<1 gram for small species of 
animals) also preclude the inclusion of ultra-precise clocks into transmitters which in 

15 turn prevents the implementation of sub-meter precision time-of-arrival location 
methods. 

Much of the existing work on radio location has been tailored to the CDMA based 
cellular telephony and GPS applications. Both these applications have the benefits 
of relatively high transmission power, coherent demodulation and stablity. Cellular 
20 telephony localization schemes are typically based on statistical measurement of 
signal strength, direction of arrival (DoA), time of arrival (ToA) or time difference of 
arrival (TDoA) as a means to position a radio source: Caffrey, J.J., Wireless Location 
in CDMA Cellular Radio Systems, (Kluwer academic Publishers: Norwell, 
Massachusetts, 2000). 

25 Although DoA methods such as beam forming can be applied to wildlife radio 
transmitter location finding, estimator precision will vary according to the number of 
antenna elements used and the relative location of the radio source from the antenna 
array. Moreover, unambiguous location finding by means of beamforming 
necessitates specific antenna geometries that could under field conditions restrict the 

30 applicability of the technique. 
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DoA methods based on orthogonal projection such as MUSIC (Schmidt, "Multiple 
Emitter Location and Signal Parameter Estimation," IEEE Transactions on Antennas 
and Propagation, vol. AP-34, No. 3, (Mar. 1986), pp. 276-280,) offer improved 
precision performance over classic beamforming. Such methods present two 
5 drawbacks for wildlife monitoring. First, the number of radio sources that can be 
monitored cannot exceed the number of antennas. Typically, wildlife studies involve 
the simultaneous monitoring of large numbers of transmitter tagged animals. Second, 
an appropriate array manifold selection must be made (Schmidt, "Multilinear Array 
Manifold Interpolation", IEEE Transactions on Signal Processing, vol. 40, No. 4, (April 
10 1992), pp. 857-866). This typically entails an empirical measurement of the antenna 
manifold. That is to say that a radio transmitter must* be used to calibrate the system 
by moving the device across a fine grid of points encompassing the search area. In 
most animal monitoring applications the empirical calibration of the array manifold 
would be logistically or economically impractical. 

15 Location finding based on ToA and TDoA techniques (including GPS) are geometry 
dependent (i.e. transmitter and receiver locations determine precision levels) and are 
precluded from use in small animal applications due to the inability to perform 
coherent demodulation and the current technological limitation of incorporating an 
ultra-stable clock source on-board a miniature (<1gram) transmitter. 

20 The measurement-based location-finding approach in Wax et al. in U.S. Pat. No. 
6,104,344 (which is not admitted to be prior art with respect to the present invention 
by its mention in this Background section) provides an alternative to past 
radiolocation finding methods. The approach calls for a library of covariance matrix 
related signatures to be empirically collected across the search array. This library is 

25 then compared using a statistical technique against measured signatures from 
mobile transmitters to determine location. This approach, like the empirical DoA 
manifold method outlined by Schmidt, has limited applicability to wildlife tracking from 
a logistic, portability and economic perspective. 

SUMMARY OF THE INVENTION 

30 This invention allows radiolocation finding with a precision on the order of V* of the 
signal carrier wavelength by exploiting a virtual interference pattern created by 
combining phase measurements made by an antenna array and a library of complex 
vectors describing the phase relationship between individual antenna and 
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hypothesized source locations. The random geometry of the antenna locations plays 
a key role in the invention. 

A randomly distributed array of antennas exhibits a virtual interference pattern with 
fine scale fractal-like structure. As the number of antennas in the array increases the 
5 fine scale structure reaches a pattern of uniqueness in relation to the location of a 
radio source. The uniqueness property of the interference pattern in relation to the 
source location is revealed with a sufficiently small spatial sampling aperture. 

In a preferred embodiment of the invention, CDMA enabled transmitters are located 
in the search area (reception zone) by first isolating the signal from one transmitter 

10 as received at multiple antennas. Separation of individual transmitter signals involves 
a form of replica code correlation that specifically preserves the phase information. A 
snapshot vector of phase information is formed for each detected transmitter and 
processed separately to form a localization estimate using either a direct search over 
the virtual interference pattern or by means of trajectory signature matching. This 

15 technique of signal separation allows orthogonal projection to be used when the 
number of transmitters in the reception zone exceeds the number of receiver 
antennas. An added benefit of the signal separation process is the effective reduction 
of noise power that results from CDMA coding gain. With coding gain, source 
location forgoes the need to explicitly time average snapshot vectors to compute a 

20 covariance matrix. 

In another embodiment of the invention, non-CDMA enabled transmitters can be 
located by building an interference pattern and performing a location search directly 
without first isolating the signal of interest. Two additional operations are required to 
support radiolocation of non-CDMA transmitters. First, prior to location estimation the 
25 number of transmitters captured by the snap shot covariance matrix must be 
identified. This operation defines the dimensionality of the covariance matrix noise 
sub-space and determines the number of transmitter peaks present in the 
interference pattern. There are a number of information-theoretic methods known in 
the art that can be used to perform this estimation. 

30 Second, non-CDMA operation requires a full rank snapshot vector covariance matrix 
to be computed. For CDMA transmitters this requirement is not required since the 
signal dimensionality of the snapshot covariance matrix will always equal one. 
Moreover, without coding gain, noise power can only be reduced by the explicit time 
averaging of snapshot vectors to compute a covariance matrix statistic. The period of 
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noise averaging will be dictated by the stationarity of aggregate radio sources within 
the reception zone of the antenna array. 

The covariance matrix, used by either orthogonal projection or trajectory signature 
matching, has a maximum rank or dimensionality bounded by the number of 
5 antennas in the array, M. Essentially; the covariance matrix cannot retain information 
on more than M transmitters. For CDMA radiolocation this bound has no impact on 
the ability to simultaneously monitor more sources than antennas. However, without 
the ability to perform signal separation, the number of transmitters cannot exceed M 
for non-CDMA radiolocation. 

10 Multi-path propagation is ubiquitous in radio. Prior art has used a number of 
methods in an attempt to either eliminate specular reflection from antenna 
measurements or by incorporating a multi-path model, typically through empirical 
measurement, into the location estimation process. By increasing the number of 
antennas in the array, the invention reduces the coherence of multi-path 

15 measurements thereby eliminating peaks in the interference pattern due to specular 
reflection. 

The invention is applicable to radiolocation for continuous wave (CW) and low duty 
cycle pulsed radio devices. 

DETAILED DESCRIPTION 

20 Figures 1 to 5 describe the invention with respect to orthogonal sub-space projection 
of an individual detection snapshot vector taken from a CDMA enabled transmitter. 
Figure 6 describes signature matching as an alternative to orthogonal sub-space 
projection. Figures 7 through 9 provide a generalization of both methods to non- 
CDMA signal localization. 

25 Figure 1 is a flow chart showing the front-end signal processing chain for m th antenna 
chain CDMA applications. 

Figure 2 is a block diagram of showing Complex Envelope OOK Demodulation, 
which appears as a single block in the flow chart of Figure 1 . 

Figure 3 is a block diagram showing the Correlation function used to decouple 
30 individual CDMA transmitter signals from one another. 
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Figure 4 is an expanded view of the Normalization block in the flow chart of Figure 1 . 

Figure 5 is a flow chart showing the location estimation based on pure phase angle of 
arrival information for a single pseudo random PN code decoupled from other CDMA 
signals. This figure assumes the front-end signal processing chain shown in Figures 
5 1 through 4 

Figure 6 is a flow chart showing an extension to the location estimation method of 
Figure 5 based on a trajectory of selected points on the interference pattern. 

Figure 7 is a flow chart generalizing the front-end signal processing chain for m th 
antenna chain applications system shown in Figure 1 to the non-CDMA case, where 
10 multiple transmitter signals cannot be directly decoupled from one another. 

Figure 8 shows the location estimation technique applicable in the non-CDMA case 
for using an averaged covariance matrix and a pre-computed estimate of the number 
of transmitters. 

Figure 9 shows an extension of the signature-matching embodiment to non-CDMA 
15 radiolocation. 

Figure 10 shows an example of an interference pattern generated for a CDMA 
encoded transmitter and a random array of 12 antennas. 

GENERAL DESCRIPTION 

The invention describes a system and method for locating transmitters to a known 
20 precision using a virtual interference pattern of expected phase difference of arrival 
signals from hypothetical transmitters transmitting from grid points within the 
boundaries of a random antenna array. 

To estimate source location the virtual inference pattern can be fully exposed by 
using orthogonal projection or by examining a trajectory of limited points across the 
25 pattern. The virtual interference pattern can be thought of in terms of weighting the 
measured signal at each antenna by a set of complex exponentials representing the 
relative phase delays from each point in the space. This gives, for every point in the 
space a sum of sinusoids of the form 

2>osfo -<Pkf 
k 
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Where O k and <p k are the relative phase delays at the /c* h antenna for the received 
signal and the hypothetical signal originating from spatial grid location co = (x 9 y), 
respectively. The sum of sinusoidal terms will have a maximum where the argument 
(6 k ~<p k ) is equal to any multiple of;r, or, equivalently where the relative distances 
5 to the antenna from the transmitter and the point co = (x,y) are related by an integer 
multiple of half the wavelength (including, of course, zero). The loci of maximum 
values (where the real and hypothetical "signals" interfere constructively) are circles. 
If the antenna location is (x k9 y k ) and the distance from the transmitter to the 

antenna is d k {co) = <>](x-x k ) 2 + (y-y k ) 2 , these circles have the equation: 

10 (x-x k ) 2 +(y-y k f =(d k +nA/2f, nel + 

When the interference pattern is fully exposed, the underlying structure of 
constructive and destructive interference reveals maxima, or points of magnification 
in which phase differences of arrival overlap constructively forming very narrow 
boundary effects or caustics. Such points correspond to potential source locations, 

15 with the global maximum corresponding to the location of the true signal. 
Alternatively, sampling a trajectory across a sufficient number of points in the 
interference pattern yields a unique signature vector. By comparing the measured 
signature vector to a stored library of hypothesized signatures a source location 
estimate can be made. Whether source location is extracted through searching the 

20 entire interference pattern or by generating a trajectory signature across a sub-set of 
points on the pattern, a direction angle is never computed; instead a phase signature 
is essentially compared with hypothetical signatures for all possible locations in the 
array area to yield the source position directly. 

In practice, the interference pattern is built not from the received signal itself but by 
25 treating the pattern as an orthogonal projection in the "noise subspace" of the 
measured signal covariance matrix. In other words, for an M element antenna array, 
we search the set of hypothetical signals looking for the peaks of a "spectrum" given 
by 



k=m+l 
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Where p k , k = m + 1, ...M denotes the noise eigenvectors of the measured signal 

covariance matrix. Here the dimension of the noise sub-space is taken to be M -m, 
such that m defines the number of independent signals encapsulated by the 
covariance matrix of antenna snapshot vectors. Linear arrays translate signal delays 
5 into angles of arrival (of plane waves), so that a{co) corresponds to a direction 
vector, and the parameter co is the hypothesized incidence angle. Alternatively, the 
invention operates from a non-linear array that translates phase delays into locations 
in space, where a{co) is the vector of relative phase delays of a hypothetical signal 

originating at a point a> = (x, j>). 

10 To see the connection with the virtual interference pattern, note that for the case of 
one sinusoidal signal the signal subspace eigenvector /?, will be identical with the 
baseband signal vector itself. If, instead of looking for zeros in the noise subspace, 
we were to look for a maximum in the signal subspace, the spatial spectrum would 
take the form 

15 PM = |^H 2 =K^'aH 2 

Where a(£) is the vector of delays corresponding to the measured signal. Each of 
these products, for a given value of the location estimate a) = (x,y), is a sum of 
sinusoids of the form 

M 

= y£ e j2*(d k (t)-d k (a>))/A 
k=l 

20 When co = % the terms of P(co) add in phase resulting in an interference pattern 
peak at the source location. 

The invention claims that a system can be designed to ensure with high probability 
that only one peak will occur for a given transmitter signal. For CDMA signals, there 
will be one peak per interference pattern. For non-CDMA signals, the number of 
25 peaks in the interference pattern will correspond to the number of transmitters. The 
practitioner positions an antenna array and computes a library of manifold vectors, 
one for each grid point. The array may be tested for uniqueness using the 
perturbations of hypothetical transmissions originating at points on the sample grid, 
or by evaluating a set of hypothetical signals on a "worst-case" alternate grid (offset 
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by V* of a sampling interval from the sample grid). The hypothesized transmitter 
signals are processed by the location algorithm to ensure that the test locations 
generate true maxima (at or near the closest points on the sample grid). If 
ambiguities are found, additional antennas can be added to the search area, and/or 
5 the spatial sampling rate may be increased, and the experiment repeated until a 
suitable differential exists between the hypothesized location and any near matches. 
The signature method may require a different number of antennas than the 
orthogonal projection method and a different threshold. 

DESCRIPTION OF PREFERRED EMBODIMENTS 

10 The invention can be described with reference to 4 preferred embodiments. 

The first embodiment is shown in Figures 1 through 5 where the localization of a 
CDMA enabled transmitter is effected using a single detection snapshot vector. The 
invention is applicable to all forms of radio CDMA transmission and applicable 
modulation schemes based on PN signal encoding, including CW communication 
15 signals and low duty cycle transmissions. In the example, an L-bit PN code family 
defines the transmit symbol library used in modulating a carrier by means of on-off- 
keying (OOK) signalling and a fixed code transmission rate. 

Derivation of a location estimate requires a number of signal processing stages. 
These processing stages are described for the antenna as follows. 

20 1 . An antenna array consisting of M antenna elements is organized in a random 

pattern, with the elements distributed so that the outer boundary of the array 
defines the search area. The invention claims that any random distribution of 
antenna elements within a pre-defined study area will generate a unique and 
precise estimate of transmitter position provided the number of antennas 

25 used is greater than a fixed minimum. Thereafter, further antennas will not 

increase the resolution of the system beyond 1/4 wavelength. However, 
increasing the number of antennas will reduce multi-path coherence, thereby 
eliminating peaks in the interference pattern due to specular reflection. 

2. Figure 1 shows the front-end processing stages of the antenna channel. 
30 Prior to digitization of the signal, the antenna output is band-limited, equalized 

for amplitude and phase imbalances and down-converted to an I.F. 
frequency. 
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Analog-to-digital conversion takes place at an intermediate frequency (I.F.) 
using a sampling rate of T s (seconds) to generate the real sampled 

signal, x m [n] (Figure 1). Nyquist sampling is assumed, so that all signal 

information is contained in a band located below the Y 2 sampling frequency 
point. A common master clock maintains sampling phase coherence across 
all antenna channels to prevent signal phase distortion. 

Generation of the complex envelope x m [n] (Figure 2) provides OOK 

demodulation while preserving signal phase information. Envelope generation 
is performed coherently for CW signals assuming sufficient signal-to-noise 
ratio (SNR), or non-coherently, for low duty cycle burst signaling and low SNR 
CW signals. Following demodulation, a sample decimation operation, rate K a , 
creates a sequence at two times the code rate. 

PN code replicas are stored as a library of sampled square waveforms. Each 
stored waveform consists of a 2xL sample sequence of +1 and -1 
representing binary code bit samples (Figure 3). The code waveforms are 
sampled at two times the code rate. 

When a transmitter is operating within the reception zone of the antenna 
array, a correlation operation (Figure 3) matches the input sample sequence 
with a stored code replica while preserving signal phase information. 
Correlation is performed for each PN code stored in the library. The auto and 
cross correlation properties of PN codes provide the means of code 
discrimination through coding gain. 

The correlation operation (Figure 3) results in a superposition of complex 
phase samples associated with an individual transmitter signal. For signal 
samples not matched to a stored code replica, superposition does not occur. 
When signals overlap in time, as would occur when multiple transmitters are 
inside the reception area of the antenna array, the inherent coding gain of the 
system enhances through superposition the phase information of one signal 
while suppressing phase information from other overlapped transmissions. 
This property of the CDMA embodiment creates virtual interference patterns 
for each transmitter. Moreover, the ability to isolate each transmitter in this 
manner enables the embodiment to simultaneously monitor more transmitters 
than there are antennas. 



8. The detector operation (Figure 4) creates a threshold level based on the time 
varying low-pass envelope of each correlator output. When the correlator 
output exceeds the threshold a detection event is registered. Other 
transmitted codes that have not been matched to a stored replica result in a 

5 correlator output that does not exceed the detector threshold. 

9. Environmental factors and related propagation mechanisms will cause the 
magnitude of a detection event to vary according to time and geography. The 
output of the detector for each antenna (Figure 4) is normalized to unity, 
removing propagation loss dependencies from entering the location 

10 estimation process. 

10. Following front-end signal processing (Figure 1), detection information from 
the antenna array is formed into a measurement snapshot vector for each PN 
code (Figure 5). The snapshot vector encapsulates the phase information of 
a single transmission event as measured across the antenna array. 

15 11. A null space is defined through an orthogonal partitioning of the snapshot 

vector into signal and noise only subspaces by means of singular value 
decomposition (Figure 5). The subspaces retain phase-difference-of-arrival 
information between antenna pairs of the received signal. 

12. Since detection is performed prior to the formulation of the snapshot vector, 
the dimensionality of the signal sub-space will be fixed to one, regardless of 
the number of transmitters within the reception range of the antenna array. 
(Except of course when there are no transmitters, in which case the 
dimensionality of the signal sub-space is zero.) This is a key element of front- 
end signal processing exploited by the invention. The invention can operate 
when the number of transmitters in the reception zone is variable or exceeds 
the number of antenna array elements. 

13. Location estimate computation is then performed using orthogonal sub-space 
projection (e.g. MUSIC). A location estimate is formed by projecting an array 
manifold vector (modeled as the snapshot vector corresponding to a 

30 hypothesized transmitter at a specific set of coordinates) onto the noise 

subspace of the measurement snapshot vector. The projection operation is 
performed across a pre-stored uniform set of grid point locations (Figure 5) 
defining the search area. The spacing of the grid must be less than or equal 

10 
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to 1/8 the carrier wavelength of the transmitted signal. In general the finer the 
grid, the smaller the minimum number of antennas required for a unique 
solution. A location estimate is defined by the largest amplitude projection 
across the set of grid points. This approach yields location estimates with a 
precision of ±% of a wavelength. 

Figure 6 describes the second preferred embodiment of the invention where the 
CDMA radio transmitter localization technique employs a trajectory of sample points 
on the virtual interference pattern. When a K » M sample point random trajectory 
is selected through the interference pattern a unique signature is generated. A 
transmitter location estimate is determined through comparing the resulting signature 
to a database of prototypes. Figure 6 describes the invention as it may be applied to 
CDMA radio transmitter localization applications using trajectory signature matching. 

1. The formation of the detection snapshot vector takes place in the same 
manner as described in Figures 2 through 4. 

2. A pre-selected K point sample trajectory (Cartesian coordinates 
P k = [x k ,y k ],k = 1,2,. K ) defines the trajectory vector V k , where the phase 

term is a function of the k* 1 location along the trajectory and the m 4 * 1 
antenna location. The shape of the sample trajectory is arbitrary. 

3. The number of sample points associated with the trajectory is selected to be 
greater than the number of antenna elements, K » M to ensure 
uniqueness of the trajectory signature. 

4. Signature samples for the I th PN code are generated by the product of K 
trajectory vectors with the PN code detection snapshot vector, 

S^V.dld^ k = 1,2,...*. 

5. A location estimate is formed by finding the nearest match of a prototype 
vector, S pt , with the measured trajectory signature, S t . The criterion for the 

best match is the minimum geometric distance between S pi and A 

prototype vector is required for each grid point on the virtual interference 
pattern. 
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Figures 7 through 9 describe the third (orthogonal projection) and fourth (signature 
matching) preferred embodiments of the invention as it may be applied to non-CDMA 
applications including those based on transponding architectures for both orthogonal 
projection and signature matching. 

1 . An antenna array consisting of M antenna elements are located in a random 
pattern distributed where the outer boundaries of the array define the search 
area of interest. The invention claims that any random distribution of antenna 
elements within a pre-defined study area will generate a unique and precise 
estimate of transmitter position. 

2. Figure 7 shows the front-end processing stages of the /77 th antenna channel. 
Prior to digitization of the signal, the antenna output is band-limited, equalized 
for amplitude and phase imbalances and down-converted to an I.F. 
frequency. 

3. Analog-to-digital conversion takes place at an intermediate frequency (I.F.) 
using a sampling rate of T s (seconds) to generate the real sampled 

signal, x m [n] (Figure 7). A common master clock maintains sampling phase 
coherence across all antenna channels to prevent signal phase distortion. 

4. Generation of the complex envelope x m [n] (Figure 2) provides demodulation 

to baseband while preserving signal phase information. Envelope generation 
is performed coherently for CW signals assuming sufficient signal-to-noise 
ratio (SNR). Or non-coherently, for low duty cycle burst signaling and low 
SNR CW signals. Following demodulation, a sample decimation operation, 
rate K a , is used to minimize computation. 

5. Environmental factors and related propagation mechanisms will cause the 
magnitude of a detection event to varying according to time and geography. 
The output of the demodulator is normalized to unity, removing propagation 
loss dependencies from entering the location estimation process. 

6. The measurement snapshot vector, d[n] , is formed from the output of each 
front-end processing channel (Figure 8 and 9). 
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7. For orthogonal projection (Figure 8) and signature matching (Figure 9), a 

covariance matrix estimate, R , is formed by averaging the outer products of 
N > M measurement snapshot vectors. 

8. For orthogonal projection (Figure 8) the null space computation requires an 
estimate of the number of transmitters captured within the snapshot 
covariance matrix. Prior art, based on information theoretic metrics 
(AIC,MDL), are used to estimate the number of transmitters and 

correspondingly, the signal subspace dimension of R . 

9. The null space of R is defined through an orthogonal partitioning of the 
matrix into signal and noise by means of singular value decomposition (Figure 
8). The orthogonal partitioning preserves phase-difference-of-arrival 
information between antenna pairs and transmitters inside the reception area 
of the antenna array. For non-CDMA signaling, the invention is restricted to 
cases where the number of transmitters is less than the number of antenna 
elements. 

10. Location estimate computation is performed using orthogonal sub-space 
projection (e.g. MUSIC). Transmitter location estimates are formed by 
projecting an array manifold (modeled as the signal component of the 
observation vector corresponding to a hypothesized transmitter at a specific 

set of coordinates x p ,y e ) t onto the noise sub-space of R . A detection 

threshold level is applied to the output of the projection operator to 
differentiate signal from noise (Figure 8). 

11. The system can use a signature matching approach instead of the procedure 
outlined in step 10 above. Signature matching (Figure 9) is performed using a 
pre-selected K point sample trajectory (Cartesian coordinates 
P k = [x k9 y k ] 9 k = 1,2 9 ...K ) which defines the trajectory vector V k , where the 

phase term is a function of the A* location along the trajectory and the 
/77 th antenna location. The shape of the sample trajectory is arbitrary. 

12. The number of sample points associated with the trajectory is selected to be 
greater than the number of antenna elements, K » M to ensure 
uniqueness of the trajectory signature. 
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13. Signature samples are generated by the product of K trajectory vectors with 
the covariance matrix R , S ki = V k RV k , k = 1,2,... K (Figure 9). 

14. A location estimate is formed by first computing the geometric distance 
between S pi and S then comparing the result with a detection threshold 

5 level (Figure 9). 

This invention is particularly useful in monitoring wildlife. A tag with a CDMA enabled 
transmitter is affixed to an animal and the animal is released into the wild. The tag 
collects data on the animal and broadcasts the data every couple of seconds using 
the pseudo random number code. When a transmitted signal is received, in addition 
10 to retrieving the data, the receiver stores phase information on the received signal 
that is ultimately used according to the methods described in the first and second 
preferred embodiments above. 

Figure 10 shows the interference pattern generated by a single 50MHz CDMA 
transmitter over an area of approximately 100 meters by 100 meters. The 
15 propagation environment was subject to random sources of specular multipath. An 
array of dipole antennas were randomly placed around the periphery of the study 
area. To achieve a unique position estimate the interference pattern was sampled at 
1/8 the wavelength or 0.75meters. 

Similarly, non-CDMA enabled tags can be affixed to targets and located according to 
20 the third and fourth preferred embodiments. Where orthogonal projection methods 
are used without a signature sample of K trajectory vectors, the number of targets 
tracked at any one time is limited to one less than the number of receivers. 

The present invention has been described with respect to low powered CDMA 
enabled tags in a wildlife tracking scenario. It will be appreciated by those of ordinary 
25 skill in the art that, in light of the present disclosure, the system and methods 
described allow for radio-location using the interferometry patterns, without departing 
from the intended scope of the invention. All such modifications are intended to be 
included within the scope of this disclosure. 



